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Abstract. Consider a system of independent random walks in the discrete torus with 
creation-annihilation of particles and possible explosion of the total number of particles 
in finite time. Rescaling space and rates for diffusion/creation/annihilation of particles, we 
obtain a stong law of large numbers for the density of particles in the supremum norm. The 
limiting object is a classical solution to the semilinear heat equation dtu — d^xU + f{u). If 
/(u) = , 1 < p < 3, we also obtain a law of large numbers for the explosion time. 



1. Introduction 

We consider nearest-neighbors symmetric independent random walks superposed with birth 
and death dynamics in the discrete torus. At rate one, each particle jumps to one of its 
neighbors with the same probability. In addition, if at a site there are r particles, at that site 
a new particle is created at rate b{r) and a particle is destroyed at rate d{r). 

We study suitable scaled versions of this process, that (as will be shown) converge almost 
surely in the L°°-norm to the solution of the semilinear parabolic problem 

ut = u^^ + f{u) gTx [o,r), 

u{x,Q) = Lp{x) > X G T. ^ ' 

where f = b — d is assumed to be smooth, ip is smooth and nonnegative, and T denotes the 
continuous one-dimensional torus. This equation has been widely studied in the literature, 
being used to model diverse processes in mechanics, physics, chemistry, technology, biology 
and many other areas. For instance, under certain conditions, it describes conduction in 
plasma, gas filtration and liquids in porous media, chemical reactions, processes of growth 
and migration of populations, etc. 

One of the most remarkable properties of this equation is the possible occurrence of sin- 
gularities due to the presence of the nonlinear source /. Even starting from regular data, for 
which there exist an existence, uniqueness and continuous dependence theory for short times, 
the solution may develop singularities in finite time. Although for linear evolution problems 
singularities may occur, they appear due to singularities in the coefficients or in the problem 
data, while in this case, singularities appear because of the nonlinear essence of the equation, 
and the time and space localization of them has to be determined through a careful analysis. 

In this problem, singularities appear in the simplest way: they are due to a fast increasing 
of the solution that leads the L°° norm to grow-up to infinity in a finite time T'max — -^max(v^)* 

The phenomenon is known as blow-up, and is interpreted as an abrupt change in the 
order of magnitude in the modeled quantity. It was successfully used to model, for instance, 
explosions in exothermic chemical reactions, population dynamics, fatigue cracking (in this 
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case explosion means that a microscopic crack changes its scale and becomes macroscopic, 
indicating a crack in the material due to fatigue). 

In view of this, it is important to understand the microscopic behavior of this kind of 
systems, and in particular their singularities (presence of them, space-time location, order, 
etc.). 

A well known condition on the nonlinear term / that assures the existence of solutions 
with blow-up is being convex, strictly positive in some interval [a, +00) and 



The most simple source verifying this conditions is f{s) = s^, p > 1. For a general 
description of the blow-up problem we refer the reader to the books |14l [T5] and the surveys 



Coming back to the Markov chain (the particle system), the scaling here is in the same 
spirit of [31 m El [in!, 'rh^ initial quantity of particles per site is also rescaled, different than 
in the hydrodynamic limit context [S]. 

As main results, we obtain almost sure convergence for the density of particles in the 
supremum norm for any compact time interval not containing Tmax- 

This result was obtained by Blount, [1] for / a polynomial with negative leading term. In 
that case solutions are known to be bounded for every t > and hence globally defined, see 



The first part of the proof consists in proving the result for birth rates b with compact 
support, where there is no blow-up and the solution is bounded. This follows essentially the 
work of Blount [4] , the main difference being that we consider any continuously differentiable 
/ instead of polynomials. In the second part of the proof, we couple a sequence of Markov 
chains as the one described above, with birth rates bj, where {bj)j>i is a sequence of smooth 
functions with compact support approximating b. This coupling allows us to extend the result 
proved in the first part to any smooth b, including those one satisfying ([2|). 

As an immediate corollary, the liminf of the explosion times of the discrete systems is no 
smaller than the blow-up time of the solution to the PDE The opposite inequality is 
much harder and we are only able to obtain it in some specific cases. Assuming that f{s) 
verifies Q and that d is a bounded or linear function, we prove that, for each A'^ > 1, the 
corresponding particle system explodes with probability one. Under the additional condition 
/(s) = b{s) — d{s) = s^, with 1 < p < 3 we also prove that the explosion time of the particle 
system converges in probability, as — >■ 00, to the blow-up time of the solution to the PDE. 

We remark here that if / is globally Lipschitz then there is no explosion with probability 
one. 

To the best of our knowledge, there is no previous work about the limit density of interacting 
particle systems exhibiting blow-up. For instance, [3llll[9l[ip] in this type of scaling, and pT] 
in the hydrodynamical limit context, have considered creation of particles, in some cases with 
unbounded limit, but all of them deal with processes defined for all times. 

It is worth to notice that since the particle system explodes in finite time, the expectation 
of the number of particles is infinity at any positive time. Hence, any method based on 
expectations is doomed to fail. This motivates the use of couplings. 
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2. Notation and results 

2.1. The particle system. Denote Tat = Z/(A^Z) the discrete torus with points. Fix two 
nonnegative smooth functions b, d : M-|- such that ci(0) = 0. Consider also a parameter 
£ G N, which will represent the number of particles per site in the initial configuration. We 
characterize the continuous time Markov chain 

{v{t))t>o = {m{t), m{t))t>o 

with state space f^A? = N^^ U {oo} by its jump rates given by 

• at rate N'^rj^, a particle jumps from k to k + 1; 

• at rate N'^rjk, a particle jumps from k to k — 1; 

• at rate £b{£~^r]k), a new particle is created at k; 

• at rate id{i~^rjk), a particle is destroyed at k. 

The transitions above are assumed for all k G Tat. Aiming not carry on the notation, we do 
not index rj[t) on N and on L Since there are no assumptions on the behavior of h at infinity, 
the waiting times of this Markov chain can be summable. If that is the case we say that the 
process explodes or blows up, and we define the state of the process as oo for times greater or 
equal than the sum of the waiting times, that we call T^^-^. More precisely, define 

:= inf{t > 0: ||r?(t)|U > iM] and T^,, := lim . 

Hence we can easily define r]{t) for t < T^^^ and we define r]{t) = oo for t > T^^^. A graphical 
construction of this process is given in Section [3l For more on explosions of Markov chains 
we refer to [12j. 

Next, we define the spatial density of particles , which is a function defined on the 
continuum torus T = [0, 1], identifying and 1 and considering a fixed orientation. 
For k G Tjv, let Xk = k/N and define 

x^{xk,t) = r^7]k{t). 

Also for Xk < X < Xk^i, we define the density of particles by linear interpolation, i.e. 

X^{x,t) = {Nx-k)X^{xk+i,t) + {k + l-Nx)X^{xk,t). 

If ry(t) = oo, we say that t)||oo = oo as well. We point out that this interpolation has 

no special meaning. Using instead a smoother interpolation or even defining X^ as a step 
function would not change the results. 

2.2. The partial differential equation. We make the following assumptions on the data 
of problem ([T]): 

• The initial datum ip is and nonnegative. 

• The source term / is and /(O) > 0. 

Under the above assumptions this equation has a unique (local) solution u, which is smooth 
in some interval (0, Tmax)- Here Tmax is the maximal existence time. If / is globally Lipschitz 
then Tmax = +oo (global existence) but it can be proved (see Theorem l6.1l in Section [6]) that if 
/ is convex and verifies ([2]) then, for positive the solution blows up in finite time, meaning 
that it is smooth in (0,Tmax)5 but 

lim ||'u(-,t)||oo = oo . 

ax 
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For general references on the study of this equation, existence, uniqueness and asymptotic 
behavior (including the blow-up case) see the books |T3l [TH [T5 | [IB]. 

The partial differential equation ([1]) satisfies a comparison principl43- Since we require 
/(O) > 0, n = is a sub-solution to this problem and hence for nonnegative initial data, 
the solution is positive. All the regularity assumptions on the data of the problem are to 
guarantee the smoothness of the solution. They can be relaxed to some extent, but since we 
are not focused on the problems arising due to the the lack of regularity of the solutions, we 
prefer not to include them to simplify the exposition. 

We are in position to state the main results of this paper. 

Theorem 2.1. Assume that 

(Al) \\X^{-,0) - ipi-)\\oo almost surely; 

(A2) for any OO, £ = 1{N) satisfies Y1n>o N^e'"^ < OO . 
Then, for any T G [0,rmax), 

lim sup t) — t) |[oo = 0, almost surely. (3) 

Remark 2.2. By Al, the parameter i represents the order of the initial quantity of particles 
per site. Condition A2 states that the growth of £ cannot be too slow in comparison with N. 
For instance, £{N) = N'^ satisfies A2 for any e > 0. 

An immediate corollary of Theorem 12.11 is the following 

Corollary 2.3. //Al and A2 hold, we have 

lim inf T^^^ > T^^x almost surely. 
N—^oo 

The left hand side in the above equation can be infinity in general. Next we find conditions 
to guarantee that that is not the case. 

Proposition 2.4. Let b ^ be convex, positive and such that l/fe(s) ds < +oo. Assume 
also that d is bounded or linear, then for every N we have 

P(Tiy,,<+oo) = l. 

Conditions on the growth of b cannot be removed as can be shown with simple examples. 
The convexity assumption is technical, but we are not able to remove it. It can be weakened 
to some extent assuming that b is convex on some interval of the form [a,-|-oo). Finally, we 
obtain 

Theorem 2.5. Assume d is bounded or linear and b{s) = s'^ + d{s) with 1 < p < 3. If Al 
and A2 hold, then 

lim T^^^ = Tjnax ; in probability. (4) 

The paper is organized as follows: in Section [3] we give a graphical construction of the 
process. Section [4] deals with a semidiscrete approximation of equation ([T]). We discretize the 
space variable and prove that the solution of the ODE obtained with this procedure converges 
to u as the mesh parameter goes to zero. In Section [5] we prove Theorem 12.11 We first prove 
this theorem for birth rates b with compact support relying on results of Blount [4j and then 
we extend these results to the general case that allows explosions by means of truncation and 
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a coupling argument. Finally, in Section [6] we study the asymptotic behavior of the particle 
system. We prove that the process explodes with probability one, Proposition 12.41 and that 
the blow-up times converges to the explosion time of the PDE ([1]), Theorem 12.51 

3. Graphical construction 

In this section we give the so called Harris graphical construction of the process: we 
construct the particle system as a deterministic function of a family of Poisson processes. 
This construction will be useful on the one hand to couple processes with different birth rates 
b and on the other hand to construct jointly a birth and death process that bounds the total 
number of particles in the system from below. 

Let (7V"^(i), M~{i), 1 > f) be a family of Poisson processes in with rate A^^. Let also 
{Mb{r,k), Md{r,k), r > 1,1 < k < N) be a family of Poisson processes in M+ with rates 
£b{i~^r), £d{£~^r) respectively. All the processes are taken independent. We construct a 
process ^(t) = {(,i{t),£,2{t) ■ ■ ■) that determines the position of each particle and the process 
rj(t) is defined as the empirical measure of (,{t), i.e. 

K(t) 

rj,{t) = ^l{m = k}. 

i=l 

The variable K(t) is the total number of particles in the system and will be defined inductively 
jointly with the construction of the process. Liitially we start with K{0) labeled particles 
1, . . . ,K{0) at positions Ci, . . . ,S,k{o)- Assume the process is defined up to time s > and 
proceed by recurrence. Start with s = and define 

For t £ [s,t) define ^(t) = ^(s) and then 

1. If T e M^{i) define Ciir) = C,{s) ± 1 and = ^j(s) if j / i. 

2. If T e Mbimis)) set ^K{s)+i = k and K{t) = K{s) + 1. 

3. If T G Md{r]k{s)), let j := min{i: i,{s) = k}. Set ^^(t) = iK(s){s) and i^(r) = K{s)-1. 

In words, if r E {M~^{i)yjN^ {i)) (where i is a particle in the system), then this particle moves 
to the right or to the left according to wether r € M^{i) or r G . If r E Mb{r]k{s)), a 

new particle is created at site k, and hence K{s) increases in one. Finally, if r E Md{r]k{s)), a 
particle is killed at site k. We kill the particle with minimum index and we assign this index 
to the particle with index K{s) (and then decrease K(s) in one), so that the alive particles 
are alway the ones with index 1, . . . , K(-). The process is then defined up to time r. Put 
s = T and iterate to define ^(t) and r]{t) up to time t = T^^^. Observe that T^^-^ < oo if and 
only if r](T^^^) := lim^^j^jv ||??(t)||oo = oo, and in this case, the sum of the waiting times is 
summable and equal to T^^^. 

Coupling processes with different birth rates. Observe that if we want to construct 
two different copies r], fj of the above process with two different birth rates h, b and we have 
6(s) = b{s) for s < M then we can use the same Poisson processes Af^{i), Af~{i), 1 > i, 
■^d{fi ■)) — 1 ^''^d Afbir, •) for r < M. In this sense we obtain that almost surely r]{t) = fj{t) 
for < t < T^, the first time that the process reaches the value M. 
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4. Convergence of a semidiscrete scheme 

We now consider deterministic spatial discretizations of ([T]), keeping continuous the time 
variable. The goal is to prove convergence of such spatial discretizations to the solution u of 
the partial differential equation ([1]). This result will be used as an intermediate step in the 
proof of the Theorem 12.11 

Throughout this section, we assume that the function / in ([1]) is globally Lipschitz. 

We define the semidiscrete approximation n^(t) = [ui{t), . . . ,nj^(t)) of the PDE ([1]) as 
the solution of the following ODE system: 

(5) 

<(0) = ip{xk), k E Tn. 

Proposition 4.1. Letu G C^'^(Tx [0,T]) be a positive solution of ^ andu^{t) the semidis- 
crete approximation given by ([5]). Then, there exists a positive constant C depending on the 
C^'^(T X [0,T]) norm of u such that, for every N large enough, 

sup max |u(xfc, t) — n^(t)| < CN~'^. 
te[o,r] '^eTAT 

As a consequence of Proposition 14.11 

limsup sup max |u^(t)| < oo . (6) 

N^oo te[0,T] ^^"^N 

We need the next lemma about solutions (and supersolutions) of the following ODE system 

{z'^ = N'^{zk+i-2zk + Zi.i) + C^{\zk\ + N-'^), /ceTiv, 
\^fc(0) = 0, k(^TN. 

Definition 4.1. We say that Z = {zi, . . . ,2jv) is a supersolution of ([7]) if 
(z'i^>N'^{zk+i-2zk + Zk-i) + C,{\zk\+N-'^), keTN, 

\zkio) > 0, keTN. 

Analogously, we say that Z_ is a subsolution if it satisfies ^ with the reverse inequalities. 

Lemma 4.2. Let Z and Z_ be a supersolution and a subsolution of ^ respectively, and let 
Z be a solution of ([7]). Then 

Z{t) > Z{t) > Z{t) . 

Proof. By an approximation procedure we restrict ourselves to consider strict inequalities in 
dS]). If that is not the case, we consider Z{t) = Z{t) + Cet with adequate C, and letting e \ 
yields the result. 

We prove that Z{t) > Z{t) arguing by contradiction. Assume that there exists a first time 
and A; E Tjv such that 'Zk{t^) = Zk{t^). Then, we would have 

> - z'kiU) > N'^{zk+i{t^) - Zk+i{t^) + - Zk^i{t^)) > 0, 

a contradiction. The inequality Z{t) > Z_{t) is handled in a similar way. □ 
Proof of Proposition For A; G T^r denote uj. = u{xk,t) and define the error function 

N 

efc := Wfc -Uk- 
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By means of Taylor's expansion, for k G T^, there exist G {xk,Xk+i) and G (xfc_i,Xfc) 
such that 

^fc+l — ~l~ '^x {^k ) "J^ ~l~ Uxx {^k ) 21N^ ~^ ^xxx {,Xk) i) gjjYS ~^ ^xxxx (cfe ) ^jjY4 

and 

1 1 1,1 

^fc— 1 ~ ^fc Ux{Xf^, ^) ^xxi^kj t) 21N'^ '^xxx{xki t) jyS ^^^^^i^ki 4] jy^ ' 

Summing the equations above and recalhng that u is the solution of JT]) gives 

^fc = ^^(""fc+l - 2-Ufe + Uk-l) + f{uk) - -^^{Uxxxx{ck) + Uxxxx{ck)). 

Writing := ^{uxxxx{ck) + i*a;a;a;a;(cfc)), we get that the error function satisfies, for k G Ttv, 

e'k = N\ek+i - 2ek + et-i) + f{u^) - f{uk) - a^N-^ . 
Since / is globally Lipschitz, there exists a positive constant C* independent of such that 
e'k < N\ek+i - 2ek + ek-i) + C,{\ek\ + N'^) , VA: G . 

Hence (ei,...,eAr) is a sub-solution of Consider the super-solution Z = {zi, . . . ,Z]\f) 
given by Zk{t) = e^'~'*^ /N'^. Notice that Z verifies By Lemma W?2\ 

ei{t) < Zi{t) < e'^^'^/N'^ for aU A; G Tat . 

Repeating the same arguments as before with — e^, we obtain 

\ek{t)\ < Zk{t) < e^^*'^ /N'^ for ah A; G T^v . 

This completes the proof. □ 

5. Limit for the density of particles. 

The following key estimate is obtained by Blount in [1]. 

Theorem 5.1 (Blount, [4J). Assume b is Lipschitz continuous with compact support. Then 
there exist constants K,a > depending on T and ip and a process Y'^ such that 

llX'^W-n^Wlloo < (l + Kte^*)f^ll^'^(0)-n^(0)||oo+ sup 



Moreover, 

P(e-^^ sup ||y^(t)||oo >e) <4A^V""'^ (9) 

tG[0,T] 

The proof in [1] considers the case where / is a polynomial with negative leading term, but 
the proof can be extended to our case with no difficulty. 

Proof of Theorem \2.1[ Assume first that b has compact support. Condition Al in Theo- 
rem 12.11 means that the bound in Q is summable in N, so Borel-Cantelli's lemma implies 

supjg[o,T] 11^ i't) ~ u'^ {'t)\\oo — > almost surely. This fact, combined with Proposition 14.11 
gives us 

lim sup \\X (t) — u{-,t)\\oo = , almost surely 

^-^°°ie[o,T] 

if b has compact support. 
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For general b, we consider 

M = sup ||n(-,t)||oo, (10) 
te[o,T] 

which is finite since we are imposing T < Tmax- 

Let 6m+i be a smooth function with compact support that coincides with b in the interval 
[0, M + 1]. Denote by X'^'^^~^^{x,t) the process with creation of particles driven by 6m+i 
instead of b. 

By making use of Harris graphical construction of Section [3l for each we can couple the 
processes X^'^^^^{x,t) and X^{x,t) in such a way that their trajectories coincide up to the 
stopping time 

K+i = inf{t>0; ||X^''-'+i(t)|U>M + l}. 

Observe that 

\\X^{t)-u{;t)\\^ < -X^'*'^+i(t)||oo + -n(.,t)|U (11) 

Denote by u^'^^^{-,t) the solution of ([T|) with / = 6jv/+i — d. Since uniqueness hold for ([1]) and 
||ii(') i) |[oo < M in [0,r], we have u^^~^^{-,t) = u{-,t), for t G [0,T]. Thus, the second term in 
(llip can be replaced by \\X^'^'^'^^{t) — u^'^'^^{t)\\oo, which goes to zero as t oo uniformly 
and almost surely. Hence there exists a finite random Nq such that, for N > Nq, 

sup ||X^'^+^(t)|U < sup ||n(-,t)||oo + | = M + ^. 

tG[0,T] tG[0,T] 

Since X^ and X^'^'^^^ are coupled, the first term on the r.h.s of vanishes for N > Nq, 
concluding the proof of the theorem. □ 

Observe that for e > 0, Theorem 12.11 implies liminf7v_>oo Tcxpi{N) > Tmax — e almost surely. 
Letting e \ yields 

liminf rcxpi(iV) > Tmax , almost surely, 

proving Corollary 12.31 

6. The explosion times 

In this section we prove the finiteness of the explosion time Tj^^x every N, Proposition 
12.41 and the convergence as — )• oo. Theorem 12.51 Both proofs rely on a coupling with a 
one-dimensional birth and death process. 

The goal is to construct (jointly with X^) a process Y = {Y{t),t > 0) that dominates 
the total number of particles ||f?(t)||oo ■= ^kVk{t) from below almost surely and for every 
time. According to the model, the rate at which a new particle (somewhere) is created is 
given by X^fceTjv °^ rates of each individual site. Since b is convex and 

nondecreasing, 

E »(f ) = » E >(f ) 2 ^Hw E "•) >- Hw) =■ ''I"!- i"i + ■ 

(12) 

Analogously, if d is bounded, the rate for annihilation of a particle is the sum of the rates in 
each site, which we bound (for t] ^ 0) hy 
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Then, we can construct a process {Y{t)) with rates q jointly with (?7(t)), in such a way that 
||?7(i)||oo > ^(i) almost surely. We need to slightly modify the construction of ry(t). The 
construction is almost the same with the only modification that for each k, we construct the 
Poisson processes 7Vrf(r, /c), r > 1 as independent thinnings of a Poisson process Md{k) with 
rate £||(i||oo and we add (independent) uniform marks to the points of the Poisson processes 
7\4(r, /c), i.e.: to each point t € Afb{r,k) we attach a random variable Ut with uniform dis- 
tribution in [0, 1]. All of them independent of all the processes. We construct Y recursively. 
Assume that the process is defined up to time s. Let 

T(s,y(s)) :=inf{t >s: t e UkeTN^fb{Vkis),k) U Mdik)} . 

For t G [s,t{s,Y{s))) define Y{t) = Y{s) and then 

1. If T G UfeeT^A4(%(s),A:) and C/.(s,y(s)) < ^iV&(M) / E^eT^ ^^(f ) , set Y{t) = 
Yis) + l. 

2. If T G Uk^Tr^^dik), set y(r) = Y{s) - 1. 

3. Otherwise set Y{t) = Y{s). 

The process is then defined up to time r. Put s = t and iterate to define Y{t) as long as 
possible. This construction guarantees that if for some time s > we have ||??(s)||oo > Y{s) 
then 

Mt)\\oo>Y{t), foralH>s. (14) 

It is straightforward to check that in addition, y is a birth and death process with rates 
q{r,r + l) = iNb{r /iN), and q{r, r — 1) = £A^|ld||oo- The next step is to proof that Y explodes 
in finite time and estimate the explosion time of Y{t). We begin with some considerations 
about birth and death processes with rates b,. and d^, respectively. For such a process, let r^. 
be the hitting time of the state r G N and denote 

the expected time to hit r + 1, starting at r. By the Markov property, 

= u ! ^ + T^f^ ' for r > . 

Or + dr br 

Notice that X^^^.^ fr is the expected time spent by the process starting from tq G N before 
reaching oo. With an inductive procedure we derive the formula 

n n—1 ^ n—1 r , n—2 ^ n—1 j , 

r=l r=0 r=Oj=0 r=0 j=r+li=r 

We invoke now that b,. = £Nb{j^), dr = ^A^|[d|[oo- Notice also that 

"^^ ds 

< oo . 



1 b{^) 

If d is bounded (or linear), then b,./dr — > as r ^ oo and hence the three terms in (jlSp are 
finite as n — )• oo. The first one behaves as j 1/6, each term of the sum in the second term can 
be bounded by Ce~^' and using these two facts, we bound the third term. We have proved 
Proposition 12.41 for d bounded. Observe that if d is not bounded but linear, equation (jl3p 
takes the form 

^ =d(|r?|) =:qi\v\,\rj\-l). 



N 
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Hence we don't need the thinnings to couple Y and r]. Both processes move jointly to the left 
always. The rest of the proof follows along the same steps. 

The proof of Theorem 12.51 is based on a more delicate analysis of (jl5p , but we first need 
some knowledge on the solutions of ([1]). That is the context of the following theorem. 

Theorem 6.1. Assume f is nonnegative and continuously differentiable. 

(1) There exists a time Tmax > (possibly infinite) such that there exists a unique maximal 
solution to ^ in [0, Tmax)- 

(2) //Tmax < +00 we have 

lim OlU-m = +00. 

max 

(3) If f is convex and positive for r > and <f{x) > tq for all x E T, then u blows up in 
finite time Tmax CLnd moreover the following estimate holds 

Tmax — / r ds . 



f{s) 



(4) // f{u) = uP with l<p<3, then 



lim ||n(t,-)||Li(T) = +00- 

i/ J- max 

Proof. Existence and uniqueness of a smooth classical (maximal) solution up to time Tmax is 
proved in [TH] . See also [3 [HI H [IS] . 

To prove (2) consider y, the solution of the ODE y' = f{y), y{0) = \\ip\\Loo(jy By compar- 
ison arguments, we get u{x, t) < y{t) for every t and hence the blow-up time Tmax = Tma.xiv') 
is bigger than the one for y. So, assume 

limsup \\u{t, •)||loo(t) < C . 

Then, there is an increasing sequence of times (tfe) such that tk Z' Tmax and such that the 
solution to ([1]) with initial data u{-,tk) is defined in some interval [0,T], where T depends 
only on C (and not on k). Then we can extend u for times t > Tmax which contradicts (1). 

For (3) consider <I>(t) := ||^t(-,i)||Li(T). Observe that v{x,t) = (p{x) verifies vt < Vxx + fiv) 
and then the solution u with initial data verifies u{x,t) > v{t,x) = (p{x) for every x G T, 
t > 0. Differentiating and using Jensen's inequality we get 

j^m = ^ fHx. t)) dx > fim) , m = Mlhj) ■ 

Thus, $(i) > z{t), the solution of the ODE z' = f{z),z{0) = ||</'||li(t)- Integrating this 
equation we obtain 

■ds>t. 



ziO) f{s) 

Let Tz be the maximal existence time for z. Taking the limit t Tj, we first observe that 
Tz < oo and z{Tz) = +oo, and next 

Tmax ^ T^ ^ / —T-: r ds . 



, f{s) 
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For (4) we use the following dichotomy proved by Velazquez [T7j. Assume that u{x,t) is a 
positive solution to ut = Uxx + for x G {—R,R) and t € (0,T) which blows up at t = T. 
Assume also that its blow-up set 

B = {x £ {—R, R) ■ limsup \\u{t, = +00} 

t/'T 

is contained in [—6, 5] for some 6 < R. Then B is isolated and for any blow-up point x £ B 
one of the following holds. 

|a; — \~ ( 8p \~ 



lim_ ■ — - u(x,rmax) 



log|x — x||y ' — 1)^ 

m 1 

lim_ \x - x|p-lu(x,rmax) = {{P " 1)C) , 



X /'X 



where C and m are positive constants with m > 4. Also Chen and Matano proved [6] that 
if the initial data cp is non-constant, the number of blow-up points of ([I]) is finite (moreover, 
it does not exceed the number of local maximum of (p) and hence we can apply the above 
dichotomy in a neighborhood of a blow-up point. Fatou's lemma and straightforward com- 
putations leads to liminit/'T \W{-,t)\\ii(j-^ > \\u{-,T)\\ii(j-^ = +00 if 1 < p < 3 for any of the 
alternatives. If ip is constant, the conclusion is immediate and holds for every p > 1. □ 

Proof of Theorem \2.5[ If T^ax = 00, Corollarv 12.31 implies the theorem. We observe that this 
only happens if / has a root at tq and ip = tq. For Tmax < 00, we assume d is bounded (if d 
is not bounded but linear, the proof is similar). From (jl5p we get 

n n—1 ^ n~l r , n—2 n— 1 j , 

E - E E E E 

r=ym r=ylN r=yiN j=0 ■' r=ym ^ j=r+li=r 

^ n—1 n—1 n—2 ^ n—1 

< ^ V - + hC V e-'- + — V V e-(^-^) (16) 

- IN ^ hir/m) ^JO^ ^ ^ IN bir/m) ^ 

r=ylN ^ ' ' r=yiN r=ylN ^ ' ' j=r+l 

f°° 1 , yiN „ f°° 1 

^ / Tr^ds + C'e-—+C" / —ds. 
Jy b{s) Jy b{s) 

Since there exists ro such that Ar/hr < e^^ for r > ro and the indexes in the sums start at 
r = ylN , the constants C, C, C" are independent of I and A'". 
Fix 5 > and choose M large enough in order to guarantee that 

■ds + C'e-f+C" / -^ds < 5. 



IM b{s) Jm Hs) 

Observe that M does not depend on A^. By Theorem 16.11 there exists a time T < Tmax such 
that 

II^(-,2-)IIli(t) >^/ + i- 
Take e = 1 and apply Theorem 12. H to get the existence of a finite random A'o G N such that, 
for N > No, \\X^iT, ■) - u{T, OIU < 1, which implies ||A^(r, OIIlHt) > M. Hence, 
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For each N and from time T on we construct the process Y, with initial data Y(T) = MiN 
jointly with ri(t). Using that the expected time for the explosion time of the process Y is 
given by 

oo 
r=Mm 

and ()16p we get that the the explosion time of Y, Ty, has expectation bounded by T + 5 and 
that for N > No, this time bounds from above the explosion time T^^^^ of rj. Hence, for 7 > 
we have 

<P(Ty > T + 7) + F{N < No) . 
Hence, the finiteness of A'^o and Markov inequality gives us 

limsupP(r^^ > T^ax + 7) < - • 
Af^oo 7 

Since 5 is arbitrary and Corollarv 12.31 implies the reversed inequality, the proof is completed. 

□ 

Acknowledgements 

We want to thank Pablo Ferrari, Milton Jara and Mariela Sued for fruitful discussions. 
PG is partially supported by UBACyT 20020090100208, ANPCyT PICT No. 2008-0315 
and CONICET PIP 2010-0142 and 2009-0613. 

References 

[1] L. Arnold and M. Theodosopulu. Deterministic limit of the stochastic model of chemical reactions with 

diffusion. Adv. in Appl. Probab., 12(2):367-379, 1980. 
[2] Catherine Bandle and Hermann Brunner. Blowup in diffusion equations; a survey. J. Comput. Appl. 

Math., 97(l-2):3-22, 1998. 

[3] Douglas Blount. Comparison of stochastic and deterministic models of a linear chemical reaction with 
diffusion. The Annals of Probability, 19(4): 1440-1462, 1991. 

[4] Douglas Blount. Law of large numbers in the supremum norm for a chemical reaction with diffusion. The 
Annals of Applied Probability, 2(1):131-141, 1992. 

[5] Xu-Yan Chen. Uniqueness of the tj-limit point of solutions of a semilinear heat equation on the circle. 
Proc. Japan Acad. Ser. A Math. Sci., 62(9):335-337, 1986. 

[6] Xu-Yan Chen and Hiroshi Matano. Convergence, asymptotic periodicity, and finite-point blow-up in one- 
dimensional semilinear heat equations. J. Differential Equations, 78(1):160-190, 1989. 

[7] Victor A. Galaktionov and Juan L. Vazquez. The problem of blow-up in nonlinear parabolic equations. 
Discrete Contin. Dyn. Syst., 8(2):399-433, 2002. Current developments in partial differential equations 
(Temuco, 1999). 

[8] Claude Kipnis and Claudio Landim. Scaling limits of interacting particle systems. Grundlehren der math- 

ematischen Wissenschaften ; 320. Springer, 1999. 
[9] Peter Kotelenez. Law of large numbers and central limit theorem for linear chemical reactions with diffu- 
sion. The Annals of Probability, 14(1):173-193, 1986. 

[10] Peter Kotelenez. High density limit theorems for nonlinear chemical reactions with diffusion. Probability 
theory and related fields, 1988. 

[11] Mustapha Mourragui. Hydrodynamic limit for a jump, birth and death process. (Limite hydrodynamique 
d'un processus de sauts, de naissances et de morts.). 1993. 

[12] J. R. Norris. Markov Chains (Cambridge Series in Statistical and Probabilistic Mathematics). Cambridge 
University Press, July 1998. 

[13] C. V. Pao. Nonlinear parabolic and elliptic equations. Plenum Press, New York, 1992. 



A PARTICLE SYSTEM WITH EXPLOSIONS; LAW OF LARGE NUMBERSFOR THE DENSITY OF PARTICLES AND THE BLOA 



[14] Pavol Quittner and Philippe Souplet. Superlinear parabolic problems. Birkhauser Advanced Texts: Easier 
Lehrbiicher. [Birkhauser Advanced Texts: Basel Textbooks]. Birkhauser Verlag, Basel, 2007. Blow-up, 
global existence and steady states. 

[15] Alexander A. Samarskii, Victor A. Galaktionov, Sergei P. Kurdyumov, and Alexander P. Mikhailov. 
Blow-up in quasilinear parabolic equations, volume 19 of de Gruyter Expositions in Mathematics. Walter 
de Gruyter & Co., Berlin, 1995. Translated from the 1987 Russian original by Michael Grinfeld and revised 
by the authors. 

[16] Juan Luis Vazquez. The porous medium, equation. Oxford Mathematical Monographs. The Clarendon 

Press Oxford University Press, Oxford, 2007. Mathematical theory. 
[17] J. J. L. Velazquez. Local behaviour near blow-up points for semilinear parabolic equations. J. Differential 

Equations, 106(2) :384-415, 1993. 
[18] Pred B. Weissler. Semilinear evolution equations in Banach spaces. J. Fund. Anal., 32(3):277-296, 1979. 

Universidade Federal da Bahia, Salvador, Brazil 
E-mail address: tertu@impa.br 

Departamento de Matematica, Fag. Cs. Exactas y Naturales, UBA and IMAS-CONICET, 
Buenos Aires, Argentina 

E-mail address: pgroismaQdm.uba.ar 



